Heat capacity of a thin membrane at very low temperature 
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We calculate the dependence of heat capacity of a free standing thin membrane on its thickness 
and temperature. A remarkable fact is that for a given temperature there exists a minimum in the 
dependence of the heat capacity on the thickness. The ratio of the heat capacity to its minimal value 
for a given temperature is a universal function of the ratio of the thickness to its value corresponding 
to the minimum. The minimal value of the heat capacitance for given temperature is proportional to 
the temperature squared. Our analysis can be used, in particular, for optimizing support membranes 
for microbolometers. 
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a. Introduction: Thin free standing membranes are 
extensively used for sensing and detecting, in particu- 
lar, for mounting of microbolometersji Thermal and heat 
transport properties of such membranes are very impor- 
tant for sensitivity of such bolometers and their time re- 
sponse. At low temperatures the wavelength of thermal 
phonons responsible for the heat capacity and conduc- 
tance can exceed the membrane thickness, b. In this 
case the vibrational modes significantly differ from those 
in bulk materials^ In particular, the lowest vibrational 
mode has a quadratic rather than linear dispersion law£ 

The contributions of these modes to the low-frequency 
density of vibrational states increase with decrease of the 
membrane thickness. As a result, the low-temperature 
heat capacity also increases. As the thickness increases, 
the heat capacity is determined by higher modes having 
an essentially linear dispersion law. Consequently, the 
heat capacity crosses over to that of a bulk material. As 
a result, the thickness dependence of the heat capacity of 
thin membranes is nonmonotonous having a sharp mini- 
mum at some optimal, temperature dependent thickness. 
This minimum is similar to the minimum in thickness de- 
pendence of the ballistic heat transfer (power radiation) 
predicted in Ref. 0. However, the concrete values of the 
optimal thickness for the heat capacity and the ballistic 
heat transfer are different. Consequently, a proper choice 
of the membrane thickness can be used for optimizing of 
the time response of microbolometers mounted on thin 
free standing membranes. The present paper is aimed at 
the theory of low temperature heat capacity of thin free 
standing membranes. 

b. Vibrational spectrum: The vibrational modes of 
a thin membrane are superpositions of bulk longitudi- 
nal and transverse modes, their relative weights being 
determined by boundary conditions - both normal and 
tangential stresses should vanish. The eigenmodes are 
classified as symmetric (SM) and antisymmetric (AM). 
Both are superpositions of the longitudinal and trans- 
verse bulk modes with wave vectors k J and k', respec- 



tively. The relations between k ; and k* are 3 - 
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for SM and AM, respectively. Here k± and fen denote 
perpendicular and parallel components of the wave vec- 
tors with respect to the membrane. Inserting the dis- 
persion laws of the bulk modes, k l = u)/ci and fe* = ui/c t 
where q and Ct are speeds of transversal and longitudinal 
sound, into Eqs. |T]) and ([2]) one obtains transcendental 
equations for the dispersion laws, w Sj „(fe|i), of different 
vibrational branches. Here the subscript a stands for the 
branch type while n stands for its number. 

In addition to AM and SM there exists a horizontal 
shear mode (HS), which is a transversal wave with both 
displacement and wave vector parallel to the plane of the 
membrane. The HS mode dispersion law is 3 - 
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where n is an integer number. SM, AM and HS modes 
are the only vibrations that can exist in a membrane. 
The modes possess a very important property: as follows 
from Eqs. (Q]), ©, and ([3]), the frequencies of all modes 
scale asi 
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Figure 1 shows the six lowest branches. Only the three 
lowest branches arc gapless, consequently only they con- 
tribute to the heat capacity of a membrane at fcgT <C 
7iu>hs,i(0). The lowest HS and SM are linear at k <C b^ 1 
while the AM dispersion law can be approximated at 
small feii as^ 
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position of this minimum depends not only on the prop- 
erties of the membrane material but also on temperature. 
It is important to underline that we calculated the heat 
capacity of a membrane per unit area. To find the spe- 
cific heat capacity c v we should divide the heat capacity 
per unit area C by the thickness b. 

For temperature T -C tmctjk B b the heat capacity of a 
membrane is dominated by the contribution of the lowest 
AM branch. So we can find a crude estimate for the posi- 
tion of the minimum, 6 m i n , by equating the heat capacity 
per unit area of a bulk sample to the contribution of the 
lowest AM branch to the heat capacity per unit area of 
a thin film. 

We find the total heat capacity per unit area of a bulk 
sample by multiplying the thickness b of the membrane 
and its specific heat capacity c Vl — 



FIG. 1: The six lowest branches of vibrations in a membrane. 
The lowest AM branch is quadratic around zero. 
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The lowest AM branch is a flexural wave and the dis- 
persion law of this branch can be obtained from the 2D 
analog of Bernoulli-Euler theory— We will reproduce the 
derivation here since we will need it for analysis of the 
applicability range of our theory. Choosing the z-axis 
perpendicular to the membrane we have£ 
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A 2 u = 0, D = 
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where u is the displacement perpendicular to the mem- 
brane, p is the density of the membrane material, E is 
the Young modulus, a is the Poisson ratio and A is the 
2D Laplace operator. Searching for a solution in the form 
u ~ e l ( k n r - wt ) we obtain the dispersion law for flexural 
waves: 
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Substituting the conventional expressions for the Young 
modulus and Poisson ratio through the sound velocities 
ci and ct and density p we arrive at the same dispersion 
law as given by Eq. ([5]). This equation is valid only in 
the long wave approximation, A> 6. 

c. Heat capacity of a free standing membrane: We 
compute the heat capacity per unit area from the general 
equation 
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The contribution of the lowest AM with the dispersion 
law (0 can be easily calculated from Eq. ||5J) as 
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From the equality c v b = Cam o wc obtain 
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For the temperature 0.1K the estimated position of the 
heat capacity minimum 6 m i n ~ 400 nm. 

A more accurate procedure is based on the exact ex- 
pression jS]). Making use of the scaling relation (|4]) one 
can cast this equation in the form 
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Here £ = bk\\ and we assume that the temperature is well 
below the Debye temperature so that we can integrate 
from to infinity. One sees that the heat capacity per 
unit area is a nonmonotonous function of the film thick- 
ness, b. Its minimum can be determined by equating the 
derivative dC/db to zero. It leads to the relation 
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where A is the area of the membrane. 

The heat capacity per unit area C of a thick mem- 
brane (b 3> fvK(3c t ) increases linearly with the thickness 
b as it should for a 3D body. But for a thin membrane 
(b < KirpCt) there is a minimum of heat capacity. The 



where z* is determined from the equation J-'(z*)z* = 
—2J r (z*). The value of z* depends on the ration c t /ci 
and for ct/ci = 1.6, z* ps 1.9. 

Numerical calculations including the 30 lowest 



branches give the value b„ 



240 for T = 0.1K. which 
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FIG. 2: The heat capacity per unit area C of a membrane, 
the response time r and the derivative of the radiation power 
with temperature p as a function of its thickness b. b mill iS 
the thickness for which the heat capacity reaches its minimal 
value. We used the values ct = 6000 m/s, ci/ct = 1.6, which 
are approximately the same as in the experiments^ on SiN^ 
membranes. 
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FIG. 3: Schematic of a free membrane with mounted detector. 



In this case in the equation of motion ^ an additional 
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term appears: 
d 2 u 

P ~dW = ~ 12(1 - a 2 ) dx 4 
the dispersion law being 



T 



d 2 u 
dx 2 



Eb2 fc 4 , l k 2 



(17) 



(18) 



is not far from the crude estimate obtained using only the 
lowest mode. Substituting (fTl)) into Eq. ([P2"]) we obtain 
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The result can be summarized in the universal form 
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where 6 m i n (T) is given by Eq. Q14p . This is shown in 
Fig.H 

d. Applicability range of the theory: There are sev- 
eral phenomena limiting the increase of the heat capacity 
with decrease of the film thickness. The obvious limit is 
that the elasticity theory is not applicable when the thick- 
ness is comparable with interatomic distance. Below we 
will discuss another limit posed by anharmonic effects. 
The amplitude of thermal vibrations increases with de- 
creasing thickness. As a result, in a very thin membrane 
the harmonic approximation used in this paper becomes 
invalid. 

To make an estimate, let us assume for simplicity that 
thickness of the membrane b <C & m i n and take into ac- 
count only the AM,0 branch. Due to the thermal vibra- 
tions, the membrane will not be perfectly flat but have a 
shape described by some function f(x,y). The resulting 
increase in the area of the membrane corresponds to some 
average tension, T, of the membrane. Let us therefore 
first study the membrane in the presence of an external 
tension, and later insert the thermal value of this tension. 



When the second term in pT|) becomes greater than the 
first one the dispersion law changes to linear and the 
heat capacity of the membrane approaches some limiting 
value. Since thermal vibration produce an average strain 
they will change the dispersion law. So for a crude esti- 
mate we will evaluate the average strain produced by the 
thermal vibrations and compare its contribution with the 
first item Eb 2 k 4 /12p(l — a 2 ). We can express variation 
of the membrane area as AA as | J A | V/| 2 dxdy. Substi- 
tuting for / thermal modes of a square membrane and 
performing thermal average one obtains an estimate for 
the typical elongation: 
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This relative elongation strongly depends on the mem- 
brane thickness. Substituting T ps E(AL/L) in Eq. fTTp 
and comparing its contribution with the first item for 
thermal wave vectors we obtain an estimate of the criti- 
cal thickness below which out theory is not applicable: 
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Numerically, this estimate gives the thickness of several 
atomic layers, so it does no pose additional limitations 
comparing to general range of applicability of the elas- 
ticity theory. 

e. Sensitivity and response time of the bolometer: 
Suppose we use a thin free standing membrane with a 
detector on it as a sensitive element of a bolometer. In 
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presence of incident flow of energy, Pi, see Fig. |3l we can 
obtain the temperature of a membrane as a function of 
time: 

Pi(t) = P 1 (T) - P 2 (T ) + C mem (T)(dT/dt) (21) 

where C mcm is the heat capacity of the membrane under 
the detector, P2(Jo) is the heat flow from surrounding 
media to the membrane under the detector, P\{T) is the 
heat flow from the membrane under the detector and To 
is the temperature of surrounding media. We assume that 
phonons are ballistic^ i. e., the mean free path of the 
phonons is much longer than the size of the detector. 
Usually the relative temperature difference is less than 
0.01, so we can expand the radiation power Pi(T) around 
To and use heat capacity as a constant. In this case we 
obtain a response time r = C/(dP/dT). In equation (jXSJl 
the part that depends only on bT can be separated 



C = T 2 C(bT) 
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The equation for power radiation^ can be similarly trans- 
formed to 



P = T 3 P(bT) . 



(23) 



From this we can obtain response time is a function only 
of bT 
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From this we immediately obtain that the minimal re- 
sponse time, r m i n , as a function of membrane thickness 
does not depend on temperature. The dependence of the 
response time r versus the membrane thickness is shown 
on the Fig. [2] The position of T m j n is close to the posi- 
tion of C m in- On Fig. [2] this difference is difficult to see 
because it is only 5%. If we measure pulses of energy 



shorter than r the rate of temperature increase depends 
only on the heat capacity of the membrane. In the op- 
posite case when pulses are longer than r the final tem- 
perature difference depends only on p = (dP/dT). The 
dependence of p on the ratio 6/6 m in is also presented on 
Fig. [2] From this picture we can see that the heat capac- 
ity C , the derivative of the radiation power with temper- 
ature p and response time r reach their minimal values 
at different thicknesses. This difference is about 10%. 
If one want to construct a detector he wants to obtain 
either the shortest response time or the greatest tem- 
perature difference. Our calculations demonstrate that 
the best sensitivity and the shortest response time can 
be obtained by choosing the supporting membrane with 
thickness & m ; n - 

/. Discussion and conclusions: We have studied the 
heat capacity of a thin membrane at low temperatures 
such that the typical wavelength of thermal phonons is of 
the order or smaller than the thickness of the membrane. 
Because of the quadratic dispersion law of the lowest vi- 
brational branch, the heat capacity per area will increase 
with decreasing thickness below a certain thickness 6 m i n . 
The thickness of minimal heat capacity is temperature 
dependent, 6 m ; n oc 1/T. The shape of the curve C{b) has 
the universal form (|15| . If we want to use the membrane 
for support of microbolometers, the reduction of the heat 
capacity is important for sensitivity of the detector for 
short pulses, and the thickness should be chosen equal to 
frmin at the operating temperature of the bolometer. 
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